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Abstract 

We investigate the motion of a tagged spin in a ferromagnetic Ising chain evolv- 
ing under Kawasaki dynamics. At equilibrium, the displacement is Gaussian, with 
a variance increasing as At 1 ! 2 . The temperature dependence of the prefactor A 
is derived exactly. At low temperature, where the static correlation length £ is 
large, the mean square displacement grows as {t/^ 2 ) 2 ^ in the coarsening regime, 
i.e., as a finite fraction of the mean square domain length. The case of totally 
asymmetric dynamics, where (+) (resp. (— )) spins move only to the right (resp. to 
the left), is also considered. In the steady state, the displacement variance grows 
as Bt 2 ^. The temperature dependence of the prefactor B is derived exactly, using 
the Kardar-Parisi-Zhang theory. At low temperature, the displacement variance 
grows as t/i 2 in the coarsening regime, again proportionally to the mean square 
domain length. 
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1 Introduction 



The purpose of this work is to investigate the random displacement of a tagged spin in 
an Ising chain evolving under Kawasaki dynamics £Q from a disordered, unmagnetized 
initial configuration. Many features of the one-dimensional kinetic Ising-Kawasaki model 
are by now well understood El HI El IE] , although its complete analytical description is 
still lacking. In Kawasaki dynamics, only pairs of opposite spins are flipped, so that the 
magnetization is locally conserved. As a consequence, the motion of an individual spin 
can be traced throughout the history of the system. 

We aim at characterizing the asymptotic subdiffusive nature of this displacement, 
both at equilibrium and in the low-temperature coarsening regime, first for symmetric, 
then for asymmetric dynamics. We use both analytical reasoning and numerical simu- 
lations. In particular, identifying (+) (resp. (— )) spins with particles (resp. holes), the 
Ising-Kawasaki chain is recognized as a one- dimensional stochastic lattice gas, where 
particles experience hard-core repulsion, as well as short-ranged attractive interactions 
originating from the Hamiltonian (jl.lj) of the spin system. This identification allows us 
to borrow concepts and methods from the realm of interacting particle systems |H1 El 
[HI1E1E]- Fi na lly ; this study gives some hindsight into the phenomenon of cage effect, 
a subject of much current interest in glassy systems [T3|. 

Consider a ferromagnetic chain of Ising spins a n — ±1 (n = 1, . . . , N), with Hamilto- 
nian 



and periodic boundary conditions. Since individual spins can be traced throughout the 
evolution, a label k = 1, . . . ,N can be attributed to each spin. The algebraic displacement 
rk(t) of spin number k at time t has an obvious definition: it is initialized at rfc(O) = 0, 
and then increased (resp. decreased) by one unit whenever spin number k moves one 
step to the right (resp. to the left). This definition of displacement wraps around the 
sample, thanks to the periodic boundary conditions. We shall be mostly interested in 
moments of the displacement, defined as averages over all the thermal histories, for a 
given tagged spin A; in a given initial configuration. This is indeed the right procedure in 
the asymmetric case. We shall successively consider symmetric (Section 3) and asymmet- 
ric (Section 4) Kawasaki dynamics, both at equilibrium and out of equilibrium. Section 5 
contains a discussion, and a summary of our main results, listed in Table El 

2 Ising-Kawasaki chain 

Let us summarize the main features of the model which will be needed hereafter. 
2.1 Symmetric dynamics 

We choose the heat-bath rule to define the stochastic dynamics, because of its peculiar 
property that the steady-state is independent of the applied electric field in the asym- 
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metric case (HI |5] (see below). For a pair of opposite spins (a n + a n+ i = 0), the move 
(a n — > —a n ,a n+ i — > — <r n +i) is realized with probability 

W{8H) = 7^77—— , (2.1) 

exp(pd)rt) + 1 

where <57i = 2(cr rt _ 1 o rt + oYi+iOn+2) = 0, ±4 is the energy difference between the con- 
figurations after and before the move. The heat-bath rule (|2.1|) obeys detailed bal- 
ance with respect to the Hamiltonian at temperature T = 1/(3, i.e., W(57i) = 
W(— 6TC) exp(—/35H). The moves corresponding to the three possible values of the dif- 
ference 87i are listed in Table ^ with the corresponding acceptance probabilities. 



type 


5H 


acc. probability 


moves 


condensation 


-4 


e 4/3 /(e 4/3 + x) 


- + -+->--+ + 

+ - + > + + -- 


diffusion 





1/2 


++-+^+-++ 
+ _<_►_ + __ 


evaporation 


+4 


l/(e 4/3 + 1) 


+ + ► + - + - 

++^- + - + 



Table 1: Types of moves in symmetric Kawasaki dynamics, and corresponding acceptance 
probabilities with heat-bath rule. 



At finite temperature, the three types of moves take place, and the system relaxes 
to thermal equilibrium. The equilibrium state is characterized by independent bond 
variables s n = <7 n a n+1 = ±1, with 

{—1 with probability II, , . 

+1 with probability 1 — II, 

where 

n 1 



e^ + l 

is the density of domain walls. The mean energy density is therefore E = —(s n ) = 
211 — 1 = — tanh j3, and spin correlations read 

(oWn) = (tanh/3) |n| = exp(-|ra|/£), 

where the equilibrium correlation length, £ = l/|lntanh/5| = l/|ln(l — 211) |, diverges as 

p 2/3 

s 2 

at low temperature. Finally, the reduced susceptibility reads 

X = E(^n)=e 2/3 . (2.3) 



At zero temperature, evaporation moves are not permitted. The system gets blocked 
into metastable configurations, consisting of domains of at least two parallel spins. Such 
metastable configurations cannot evolve under the sole effect of condensation and diffu- 
sion moves. A study of the statistics of the metastable configurations reached by zero- 
temperature dynamics can be found in In particular the blocking time is rather 
short, increasing as (IniV) 3 for a system of N spins. 

At low but non-zero temperature, evaporation moves take place with the small prob- 
ability 

e = - l/(4£ 2 ) « 1. (2.4) 

Starting from a disordered initial configuration, the system is first attracted into a 
metastable configuration, from where it escapes on the slow time scale et ~ t/£ 2 . For 
£/£ 2 3> 1, the system enters a self-similar coarsening regime, where the typical domain 
length grows as L ~ (t/^ 2 ) 1 ^ 3 El E] [see ()H.14j) ]. Finally, the domain length saturates 
to its equilibrium value 

L eq = — = - = e 2 ? + 1 « 2f . 

q i + E n s 

The corresponding equilibration time scales as 

r cq ~ e- (2.5) 

2.2 Asymmetric dynamics 

Assume that (+) spins (resp. (— ) spins) are positively (resp. negatively) charged. In the 
presence of an electric field, they move preferentially to the right (resp. to the left). For 
the sake of simplicity, consider the limit of totally asymmetric dynamics [HE], where only 

H > — h moves are allowed. These moves are listed in Tabled with the corresponding 

acceptance probabilities. 



type 


5H 


acc. probability 


move 


condensation 


-4 


e 4 7(e 4 ^ + 1) 


- + -+->--+ + 


conduction 





1/2 


++-+^+-++ 


evaporation 


+4 


l/(e 4 ^ + 1) 


+ + ► + - + - 



Table 2: Types of moves in totally asymmetric Kawasaki dynamics, and corresponding accep- 
tance probabilities with heat-bath rule. 

Under asymmetric dynamics at finite temperature, the driven system reaches a non- 
equilibrium steady state, where (+) spins (resp. (— ) spins) are advected with a constant 
drift velocity V (resp. — V), to be derived below [see 1)4.40 ] . All the static quantities, such 
as equal-time correlations, coincide with their equilibrium values E] for a whole class 
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of stochastic dynamical rules, including the heat-bath rule (12. lj) . This does not hold with 
the Metropolis rule. 

At zero temperature, the system again gets blocked into metastable configurations, 
consisting of domains of at least two parallel spins. The statistics of metastable configu- 
rations is different from that obtained in the symmetric case. This difference is explicitly 
apparent in the special case of restricted zero-temperature dynamics (descent dynamics), 
where only condensation moves are allowed. In this situation, analytical results can be 
obtained, both for the symmetric case |T0J E] and the asymmetric one ^H] • 

At low temperature, the system exhibits self-similar domain growth for t/£ 2 3> 1, with 
a typical domain length increasing as L ~ (t/^ 2 ) 1 ^ 2 |H E] [see ()4.8|) ]. The equilibration 
time scales as 

r cq ~ e 4 - (2.6) 



3 Symmetric dynamics 

3.1 Mean square displacement at equilibrium 

Assume that the system is at thermal equilibrium, evolving under symmetric Kawasaki 
dynamics. 

Consider first the case of infinite temperature. Then all the allowed moves have equal 
acceptance probabilities 1/2. Particles (i.e., (+) spins) hop along the chain, with the 
hard-core constraint that sites are occupied by at most one particle. The motion of holes 
(i.e., (— ) spins) follows by complementarity. The process thus defined is known as the 
symmetric exclusion process (SEP) [TUl HU H2]- The particle density p is fixed by 
the magnetization of the initial configuration: M = (a n ) = 1 — 2p = 0, hence p = 1/2. 
The displacement of individual particles is asymptotically Gaussian, with a variance 
increasing as t 1 ! 2 \T§\ I7j. 

We now show that the same holds for the system at equilibrium at any finite temper- 
ature, i.e., we predict a subdiffusive displacement of tagged particles, with a symmetric 
Gaussian dispersion profile, whose variance scales as 

{r 2 )^At 1 ' 2 . (3.1) 

The temperature dependence of the amplitude A is given in (j3.11|) . 

In order to do so, we use the Edwards- Wilkinson (EW) continuum theory [20J. 
The starting point of this approach consists in mapping the problem onto an interface 
model [2U 1221 123] j by defining the height variable 

n 

hn= a ™- ( 3 ' 2 ) 

m=-N 

The reference point n = —N does not actually matter, since only height differences will 
be considered hereafter. One has e.g. 

{{hm-K) 2 ) « \m-n\x (3.3) 
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for large separations (|m — n\ ^> £) at equilibrium, where the compressibility x is j us t 
the reduced susceptibility evaluated in (|2.3j) . Height and displacement are related by the 
following approximate relationship [2*T] |221 122] 

h n {t + r) - h n (t) = -2N n {t + r, *) « -{r{t + r) - r(t)) n . (3.4) 

In this expression, N n (t+r, t) is the net number of particles ((+) spins) which have crossed 
the bond (n, n + 1) between times t and t + r, particles moving to the right (resp. to the 
left) being counted positively (resp. negatively), while {r{t+r)—r(t)) n is the displacement 
between times t and t + r of particles in the neighborhood of site n. Particles indeed move 
together almost rigidly, as the system has a finite static compressibility \ [see (|3.3j) ]. 
In the continuum limit, the height h(x,t) »s — r(x, t) obeys the EW equation fU\ 122] 

dh ^d 2 h , . .„ 
- = D — + ,(,,,), (3.5) 

where r](x,t) is a Gaussian white noise, such that {rj(x,t)r](x' ,t')) = 2A 5(x — x')5(t — t'). 
In Fourier space, ()3.5|) yields 

%,t) = K(g, 0) e- Dg2 ' + /' dsfj(q,s)e- Dq2(t - s) . 

Jo 

We thus obtain the following expression 

fOO rlfl 1 - P - D <} 2t / T \ 1/2 

<(/>(*, t + r) - ft(i, t ))»> 2A ^ = 2A (_) 

for the height (or displacement) correlation function at equilibrium (t — > oo). The noise 
intensity A can be evaluated from the height correlation at equilibrium, 



{(h(x,t) -/i(0,t)) 2 > ^2A 



oo 



dg 1 — cos gx Ax 



2tt L>g 2 D 
Comparing with ()3.3)1 yields the identification A = Dx, and finally 

(£] 1/2 



A = 2 X • (3.6) 



We now estimate the magnetization diffusion coefficient D. A hierarchy of non-linear 
kinetic equations for spin correlation functions can be written by summing the contribu- 
tions of the moves listed in Table ^ For the magnetization, we obtain 

-j^ = S n -\ + SVi+i ~ 2S n , (3-7) 

with 

1 e 4/3 - 1 

$ri = ^ ( a n) + 4 ^4/3 _|_ (^n-l^n^+l - °n-l ~ On ~ CT n+ i). (3.8) 



6 



The nonlinear equation ()3.7|) cannot be solved in full generality, except at infinite temper- 
ature [21], where it becomes linear. However, since (J3.7|) already involves a second-order 
difference, it is sufficient for the present purpose to evaluate the uniform value S of S n 
in the finite-temperature equilibrium state of the Ising chain with a small uniform mag- 
netization M = (a n ). Indeed, f!3.7p implies 

D = lim (3.9) 

The explicit calculation of 5* can be performed by means of the transfer-matrix formalism. 
Introducing a uniform reduced magnetic field H, the transfer matrix reads 



r 



' e /3+H e -P-H 
e -(3+H e (3-H 



The leading eigenvalue of T and the corresponding eigenvectors read 

A = e^ cosh H + e~ p (l + ^ sinh 2 H) l/ \ 

(L\ = (A-e^,e~^), \R) = (/I^) . 
We have in particular 



where the spin operator reads 



hence 



1/2 



&*smhH H M+(M 2 + e 4 ^(l-M 2 

(l + e^sinh 2 ^) 1 ^' ~ e^l - M 2 ) 1 ^ ' 1 ; 

The equilibrium value of any spin correlation can be expressed in terms of matrix elements 
involving the eigenvectors (L\ and \R), and the spin operator S. Using the second 
expression of (|3.10|) . results can be written in terms of M alone. We have thus 

_ (L\STSTS\R) 

\&n-l&n&n+l) — 



A 2 (L\R) 

„xl/2 .. _ _i 

M 



8 (M 2 + e 4 ^(l - M 2 )) 1/2 + e 4 ^(e 4 ^ - 6)(1 — M 2 ) — 3 — 5M 2 



( e 4/3 _ ^2^ _ M 2) 



M (M->0). 



( e2/ 3 + 1)2 

Inserting the latter estimate into ()3.8j) . ()3.9j) results in 



Finally, using (|3.6|) . we obtain the temperature dependence of the prefactor A: 



( 



8eU 



) 



1/2 



( 



) 



1/2 



A = x 



7r(e 2 / 3 + l)(e 4 / 3 + l) 



(3.11) 



7T 



In the infinite-temperature case, we thus recover A SEP = (2/ir) 
In the low-temperature regime, the amplitude A vanishes as 



i 2 nam. 



A^(8/7r) 1 / 2 e-^2(7re)- 1/2 , 



(3.12) 



hence the low-temperature scaling form (r 2 ) ~ (t/£) 1//2 . The latter result is compatible 
with a smooth crossover, for times t ~ r oq ~ £ 5 [see f)2.5p ]. between the behavior (|3.16|) 
in the coarsening regime and the behavior f|3. 1|) in the equilibrium regime. 

In order to test the accuracy of Monte-Carlo numerical simulations, which will be 
used more extensively hereafter, we have checked the prediction (|3.11j) against values 
of A obtained by extrapolating numerical results for (r 2 ) at equilibrium, for various tem- 
peratures. Figure ^ shows a plot of A against e -13 , both for an unmagnetized equilibrated 
initial condition (initial temperature = final one), built recursively using ()2.2|) . and for 
an unmagnetized disordered initial condition (infinite initial temperature). In the equi- 
librated case, the measured values of the amplitude are in perfectly good agreement 
with (J3.ll)) . In the disordered case, A is found to be systematically smaller than the 
predicted value ()3.11|) . The amplitude A actually depends continuously on the initial 
condition, as can be observed by preparing the system at equilibrium, at a temperature 
different from the chosen final temperature. This dependence on the initial configuration 
is somewhat surprising, since the amplitude A characterizes the long-time behavior (|3.1j) 
of the displacement of a tagged particle. Let us mention that data points for lower 
temperatures are more and more difficult to obtain, because the equilibrium regime is 
reached for t ~ r eq ~ £ 5 ~ e 10/3 . 

3.2 Mean square displacement in the coarsening regime 

Consider the displacement of a tagged particle after a quench from an infinite-temperature 
disordered initial configuration to a low temperature such that £ ^> 1. 

The system is first attracted into a metastable configuration, from where it escapes 
on the slow time scale et ~ t/C 2 (see Section 2.1). Correspondingly (r 2 ) reaches first a 
plateau value, from which it escapes on the same time scale. For t/C 2 ^> 1, the system 
enters a self-similar coarsening regime, and (r 2 ) is observed to grow as {t/C 2 ) 2 ^- Finally, 
as the domain length saturates to its equilibrium value, for t ~ r eq ~ £ 5 , (r 2 ) crossovers 
to its equilibrium behavior described in Section 3.1. 

We now provide an explanation for the observed anomalous diffusion for t/C, 2 ^> 1. 
We use the picture of Cordery et al. and Cornell et al. .3] for the low-temperature 
coarsening regime (see for further developments). Consider a (+) spin evaporating 
from the right boundary of a (+) domain, of length £i, say. It then undergoes a free 
random walk inside the right neighboring (— ) domain, of length £, until it condensates, 
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1.0 




exp(-0) 

Figure 1: Plot of the amplitude A [see (|3.1|) ] of the spin displacement for symmetric equilibrium 
dynamics, against . Full symbols: measured amplitude for an equilibrated initial condition 
(initial temperature = final one). Full line: prediction (|3.11|) of EW theory. Open symbols: 
measured amplitude for a disordered initial condition (infinite initial temperature). Dashed 
line: polynomial fit to the latter data, meant as a guide to the eye. 



either (i) to the (+) domain from where it had evaporated, or (ii) to the next (+) 
domain to the right, of length £2, say. Event (i) amounts to no domain displacement at 
all, while event (ii) amounts to moving the (— ) domain by one lattice spacing to the left, 
hence £\ — *■ £\ — 1, and £2 — > £% + 1. The probability of event (ii), knowing the spin has 
evaporated, is exactly l/£. This is a well-known result in the gambler's ruin problem |2BJ. 

In other words, integrating over the fast processes (diffusion and condensation), one 
is led to an effective description of the coarsening process in terms of diffusion and 
annihilation of domains, where each domain performs a random walk independently of 
the others, with a diffusion constant 

*W «\«±- t (3-13) 

for a domain of length £. When a domain length shrinks down to zero, the two neighboring 
domains coalesce. Since a domain of length L evolves on a time scale of order t ~ 
L 2 /D(L) ~ £, 2 L 3 , this picture predicts the power-law growth [21131^] 

L^a L (t/e) 1/3 , (3.14) 

where the amplitude ai is some non-universal number of order unity. Note that the range 
of validity of the coarsening regime (1 < L C ^, or ^ 2 < f C r cq ~ £ 5 ) is larger and 



9 



larger at lower temperatures. In this regime, the distribution of domain lengths has the 
scaling form 

PW-jf(i). (3-15) 

Coming back to the displacement of a tagged particle, the scenario of domain diffusion 
and annihilation described above demonstrates that a tagged particle (i.e., a (+) spin) 
exhibits two kinds of motion in the coarsening regime: 

(a) Most of the time, the tagged particle remains inside a (+) domain. This do- 
main, and the tagged particle contained in it, perform a random walk, with an 
instantaneous diffusion constant given by (JHHEJ). These moves generate a Gaussian 
displacement, such that 




The scaling laws (ETToT) and lETTl) imply (l/£) ~ l/L ~ (i/£ 2 )~ 1/3 . As a con- 
sequence, the number of moves of type (a) and the corresponding mean square 
displacement up to time t both scale as iV( a ) ~ (r 2 ) ~ L 2 ~ (£/£ 2 ) 2 ^ 3 - 

(b) From time to time, the tagged particle evaporates from the boundary of the (+) 
domain to which it belongs, and condensates to a neighboring domain. The dis- 
placement of the particle in such a move is the length i of the (— ) domain which 
is crossed. It is therefore distributed according to the law (J3.15j) . Furthermore, for 
a given tagged particle, successive such displacements alternate in sign, because 
the particle can only bounce between two neighboring (+) domains. The num- 
ber of such moves scales as Nn^ ~ L ~ (t / 1//3 , and the corresponding random 
displacement scales as (r 2 ) ~ L 2 ~ (£/£ 2 ) 2 ^ 3 - 

To sum up, moves of either type bring similar contributions to (r 2 ). We thus predict 
the power law 

(r 2 ) « a 2 r (t/e) 2/3 , (3-16) 

where the amplitude a r is of order unity. As a consequence, the dimensionless ratio 
(r 2 )/L 2 approaches a non-trivial universal limit value 

Q= Hm !Q = (a r /a L ) 2 (3.17) 

in the late stages of the coarsening regime. The sum of moves of type (a) is asymptotically 
Gaussian. The distribution of the sum of moves of type (b) is not known a priori, although 
it might be expected to be Gaussian as well. 

We have checked the above predictions, and tested the Gaussian character of the 
distribution of the displacement, by means of numerical simulations. Figures 121 and H3 
show plots of the mean domain length L, obtained by measuring the energy density 
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i , , , , , , , , , 

2 4 6 8 10 
(t/ai/3 

Figure 2: Plot of the mean domain length L in the coarsening regime with symmetric dynamics, 
against (t/l; 2 ) 1 / 3 . Symbols: data for several values of £. Full straight line: least-square fit of 
the data (first three points excluded), confirming (|H.14jl . The slope yields ol ~ 1.34. 

E = — 1 + 2/L, and of the root-mean-square displacement (r 2 ) 1 / 2 , in the coarsening 
regime, against the scaling time variable (t/£ 2 ) 1//3 , for several values of £ 3> 1. The data 
are averaged over 1000 independent systems of 2000 spins, for times up to t = 1000 £ 2 . A 
local smoothing procedure has been applied to the raw data. The observed data collapse 
and linear behavior confirm the power laws (J3.14j) . (jH.lfij) . while least-square fits yield 

a L « 1.34, a P «0.71, Q » 0.28. 

Figure |U shows a plot of the kurtosis /C = (r 4 ) / (r 2 ) 2 , against (t/£ 2 ) -2 ^ 3 - The data 
smoothly extrapolate toward the value /C = 3 of the Gaussian law, with a correction pro- 
portional to (t/£ 2 ) -2 / 3 ~ 1/L 2 . This observation strongly suggests that the asymptotic 
law of the displacement is Gaussian. 

4 Asymmetric dynamics 

4.1 Displacement variance in the steady state 

The system is assumed to be in the steady state reached under totally asymmetric 
Kawasaki dynamics at finite temperature. 

Consider first infinite-temperature dynamics. Particles, that is (+) spins, hop to 
the right, with a hard-core constraint. This process is known as the totally asymmet- 
ric exclusion process (ASEP) [TUl HU C2]- The displacement of individual particles 
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(t/( 2 ) 1/3 

Figure 3: Plot of the root-mean-square displacement (r 2 ) 1 / 2 in the coarsening regime with 
symmetric dynamics, against {t/^ 2 ) 1 ^. Symbols: data for several values of £. Full straight line: 
least-square fit of the data (first three points excluded), confirming (|3.16j) . The slope yields 
a r « 0.71. 



is asymptotically described j2U 122] by the nonlinear Kardar-Parisi-Zhang (KPZ) the- 
ory [21]. This theory predicts that, for a given tagged spin in a given initial configura- 
tion, the displacement at time t consists of a drift with uniform velocity, and a thermal 
random displacement, increasing typically as t 1//3 . This anomalous, subdiffusive width 
is masked by usual diffusive fluctuations in t 1 ^ 2 whenever the displacement is averaged 
either over different spins or over different initial configurations. These peculiar features 
have been underlined in |2*H1 |2H] ■ 

At finite temperature, we expect the same to hold at long times, i.e., we predict that 
the motion of a tagged particle consists of a deterministic drift and a subdiffusive thermal 
motion, according to 

(r)«Vt, (r 2 ) c « Bt 2/3 . (4.1) 

The temperature dependence of the velocity V and of the amplitude B are now 
derived, using the KPZ theory. The most general KPZ equation reads 

dh dh 8 2 h A (dh\ 2 . , 

Unlike (|3.5)1 . this equation cannot be solved explicitly. However, all the quantities of 
interest can be predicted j2HI22] from the sole knowledge of the steady-state current J(M) 
for a non-zero magnetization M — 1 — 2p, where p is the particle density. This current 
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(t/f2)"2/3 

Figure 4: Plot of the kurtosis K, = (r 4 )/(r 2 ) 2 in the coarsening regime with symmetric dy- 
namics, against (i/£; 2 ) _2//3 • Symbols: data for several values of £. Full line: polynomial fit 
demonstrating a smooth convergence to the limit value /C = 3 of the Gaussian law. 



can be evaluated by summing the contributions of the moves listed in Table 121 The 
corresponding steady-state probabilities can again be derived by means of the transfer- 
matrix formalism. Let V± = (1 ± S)/2 be the projectors onto the spin states ±. The 
first move involves the probability of observing the pattern — I h: 

{L\V_TV + TV_TV + \R) 
P ~ + ~ + ~ A 3 {L\R) 

(3 + M 2 + e 4/3 (l - M 2 )) (M 2 + e 4/3 (l - M 2 )) 1/2 - (1 + 3M 2 ) - 3e 4 ^(l - M 2 ) 
= 2(e 4 ^ - 1) 3 (1 - M 2 ) ' 

Taking into account all the allowed moves, we obtain 

e 4/3 1 1 

2M 2 + e 4/3 (l - M 2 ) - (1 + M 2 ) (M 2 + e 4/3 (l - M 2 ) ^ " 
= 2(e 8 ^-l)(l-M 2 ) 

This expression shows that J(M) = J + J 2 M 2 + - ■ ■ is an even function of M, as expected 
from the spin symmetry of the dynamics, i.e., the invariance of the dynamical rules under 
the simultaneous nipping of all the spins. 

As a consequence of ()3.2j) and 1)3.40 . we are led to identify dh/dt « —2 dN/dt « —2 J 
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and dh/dx « M in the continuum limit. Therefore, in (|4.2|) we must have 

/x = -2 J , z/ = 0, A = -4J 2 . 

Furthermore, the zero-magnetization velocity V is such that Jo = pV with p = 1/2, 
hence V — 2 J = — /i. Equation ()4.3|) yields the explicit results 

e 2/3 o 2/3 _ i 

V = = yen A = — (4 4) 

{e 2 ? + l)(e^ + 1) * ' e 2 / 3 (e 2 / 3 + l)(e 4 / 3 + l) 1 ' 

According to the KPZ scaling theory [^TJ |22 EDI, for a given tagged spin in a given 
initial configuration, the displacement fluctuation is asymptotically distributed accord- 
ing to 

/X 2 AA 1/3 / e 2 ^(3e 2 ^-l) 
r {r '~\ 2 ) { 2 ( e ^ + l)(e^ + l) 

The reduced variable X has a universal distribution, related to the Tracy- Widom law of 
the largest eigenvalue of a complex Hermitian random matrix. This distribution, which 
attracted much interest in recent years (31], can be expressed in terms of a solution of a 
Painleve transcendental equation. 

As a consequence, the reduced central moments 

((r-(r>n = ( (X-Wn 

have non-trivial universal values. The variance (X 2 ) c = 0.63805, and the reduced mo- 
ments K 3 = 0.2935, K4 = 0.1652, are known numerically to a high accuracy [3T]. Finally, 
the amplitude B entering (|4.1|) reads 



2{e 2 P + l){e^ + l) 



0.63805 



In the low-temperature regime, we thus obtain 

V « e~ 4/3 « l/(4£ 2 ), B « (3/2) 2 / 3 (A 2 ). e" 4/3/3 « (3/4) 2 / 3 (X 2 ) 6; T 2/3 , (4-7) 

0.83608 0.52670 

hence (r) ~ t/£ 2 , (?" 2 )c ~ (^/£) 2 ^ 3 - These results are compatible with the existence of 
a smooth crossover between the coarsening and steady-state regimes, for times of order 

i~r cq ~e 4 [see (EH)]. 

Figure El shows a plot of the amplitude B obtained by numerical simulations, against 
e -4/3/3 J3 a t a f or lower temperatures are again very costly. The non-monotonic depen- 
dence of B on temperature is accurately described by the analytical prediction (|4.6|) 
from the KPZ theory. The data for the drift velocity V are even more accurate, and 
indistinguishable from the prediction (|4.4j) . At variance with the case of symmetric dy- 
namics, the values of V and B do not depend on the initial condition. The measured 
values indeed coincide with the theoretical predictions, up to statistical errors, both for 
an equilibrated and a disordered initial condition. 
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exp(-40/3) 

Figure 5: Plot of the amplitude B [see ([4.1)1 ] of the spin displacement in the steady-state 
regime of asymmetric dynamics, against e -4 ^/ 3 . Symbols: data. Full line: prediction (j4.fi)) of 
KPZ theory. 

4.2 Displacement variance in the coarsening regime 

We first summarize the known features of the low-temperature domain-growth process 
in the presence of a bias jlj E]. As in Section 3.2, the relevant time scale is et ~ t/l; 2 , 
since evaporation moves take place with the small probability e (cf. 1)2.4)1 ). 

A (+) spin evaporating from the right boundary of a (+) domain propagates ballisti- 
cally through the right neighboring (— ) domain, and condensates to the left boundary of 
the next (+) domain to the right. The converse situation of a (+) spin evaporating from 
the left boundary of a (+) domain is not permitted. Integrating over the fast processes 
(condensation and ballistic propagation) yields the following description of the dynamics 
in terms of domain shifting and annihilation [HE]- Each (+) (resp. (— )) domain per- 
forms a biased walk to the right (resp. to the left), with a hopping rate e, irrespective 
of its length. When a domain length shrinks down to zero, the two neighboring domains 
coalesce. As a consequence, the whole pattern of (+) (resp. (— )) domains is advected 
with a mean velocity V (resp. — V). As long as the effective description in terms of 
domain dynamics holds, the mean velocity is given by V ~ e, just as in the steady-state 
regime [see (14.7)) ] . The fluctuations around the mean ballistic motion of a given domain 
are diffusive, with a diffusion constant D ps e. 

In the self-similar coarsening regime (t/£ 2 S> 1), the distribution of domain lengths 
is still of the form 1)3.15)) . The behavior of the mean domain length L is now driven by 
the diffusive fluctuations. Since a domain of length L evolves on a time scale et ~ L 2 , 
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the following power-law growth holds [3J E] 

L « &L(t/e 2 ) 1/2 , 



(4.8) 



where the amplitude bi is some non-universal number of order unity. 

In the coarsening regime, a tagged particle therefore exhibits two kinds of motion: 

(a) Most of the time, the tagged particle remains inside a (+) domain. This domain, 
and the tagged particle contained in it, perform a biased random walk, characterized 
by a mean velocity and a diffusion constant V ~ D ~ e ~ l/(4£ 2 ). 

(b) From time to time, the tagged particle evaporates from the right end of the (+) 
domain to which it belongs, and condensates to the left end of the neighboring (+) 
domain to the right. The displacement of the particle in such a move is the length i 
of the (— ) domain which is crossed. It is therefore of order L. 

As in the symmetric case, moves of either types bring similar contributions to (r 2 ). 
We thus expect, for a given particle and a given initial configuration, a dispersion of the 
form 

(r 2 ) c ^b 2 r t/e, (4.9) 

for t ^> £ 2 , where the amplitude b r is of order unity. As a consequence, the dimensionless 
ratio (r 2 ) c /L 2 has a non-trivial universal limit value 

R= Jim ( -^ = (b r /b L r (4.10) 

in the late stages of the coarsening regime. As the displacement is biased, all the higher 
central moments K n ()4.5|) are expected to be non-trivial for finite times. If, however, 
the dispersion is asymptotically Gaussian, as it is in the symmetric case, the cumulants 
associated with the moments K n should fall off to zero in the t/£ 2 — > oo limit. 

The following results of numerical simulations corroborate the above predictions, and 
suggest that the distribution of the displacement is again asymptotically Gaussian. Fig- 
ures HI and [7| show plots of the mean domain length L and of the displacement width 
(r 2 )^ 2 , against the scaling time variable (t/^ 2 ) 1 / 2 . The statistics is similar to that of 
Figures El and H3 The observed data collapse and linear behavior confirm the power 
laws ([4.8)1 . (|4.9j) . The data for the mean domain length exhibit appreciable corrections 
to scaling, which were hardly visible in the case of symmetric dynamics. Least-square 
fits yield 

b L « 0.70, b r w 0.42, R » 0.36. 

Figure |H1 shows a plot of the reduced third central moment K$ [see (I4.5|) ] against 
(t/^ 2 )" 1 . The observed fall-off to zero strongly suggests that the dispersion is asymptoti- 
cally Gaussian in late stages of the coarsening regime. We have checked that — 3 also 
decays to zero, albeit with larger statistical errors. 
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Figure 6: Plot of the mean domain length L in the coarsening regime with asymmetric dynam- 
ics, against (t/S, 2 ) 1 ^ 2 . Symbols: data for several values of £. Dashed lines: guides to the eye, 
emphasizing the corrections to scaling. Full straight line: least-square fit of the data for £ = 30 
(first point excluded), confirming (|4.8jl . The slope yields bi w 0.70. 

5 Discussion 

As shown in this work, the motion of a tagged spin in a ferromagnetic Ising chain evolving 
under Kawasaki dynamics is akin to anomalous diffusion. The variance of the displace- 
ment, (r 2 ) c = (r 2 ) — (r) 2 , grows algebraically in time, with an exponent depending on 
the regime considered. 

In the steady state at finite temperature, the present system is in the universality class 
of well-known interacting particle systems, the SEP in the symmetric case and the ASEP 
in the asymmetric case, with subdiffusive behaviors (r 2 ) « At 1 / 2 and (r 2 ) c ~ Bt 2 ^ 3 , 
respectively. Exploiting the relationship with the corresponding continuum theories, 
EW in the symmetric case and KPZ in the asymmetric case, the explicit temperature 
dependence of the amplitudes A and B are derived exactly [see (|3.11j) and ()4.6|) ]. 

Diffusive particle systems with interactions beyond the usual hard-core repulsion have 
been the subject of some recent works. The model for polymer crystallization introduced 
in [32] is equivalent to asymmetric Kawasaki dynamics for the antiferromagnetic Ising 
chain. In this model, the coupling constant J is negative, so that (3 is replaced by —f3\J\. 
The prefactor B vanishes at a finite critical temperature such that (3\J\ = (ln3)/2, as 
noticed earlier |30]. Below the critical temperature, the current J(M) (|4.3|) exhibits a 
minimum in the unmagnetized state (M = 0), and two symmetric maxima for non-zero 
values of M such that exp(— 4/3 1 J\) = 1 — 8/(3 — M 2 ) 2 . This shape of the density-current 
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Figure 7: Plot of the displacement width (r 2 )^ 2 in the coarsening regime with asymmetric 
dynamics, against (i/£ 2 ) 1//2 . Symbols: data for several values of £. Full straight line: least- 
square fit confirming (|4.9jl . The slope yields b r ~ 0.42. 

characteristic curve is responsible for the existence of a minimal current phase in a system 
with open boundaries jHSj- The prototype model studied in |33j is also equivalent to an 
antiferromagnetic Ising chain, for a special value of the magnetic field. 

In the nonequilibrium low-temperature coarsening regime, domain growth takes place 
at the slow time scale i/£ 2 . The effective description of the dynamics on this scale is given 
in terms of domain diffusion and annihilation L3J in the symmetric case, and of domain 
shifting and annihilation |D |5] in the asymmetric case. The motion of a tagged particle 
is very similar in both cases, the variance of the displacement scaling as the mean square 
domain length. Numerical simulations yield the values Q ~ 0.28, R ~ 0.36 for the 
corresponding universal ratios [see ()3.17|) and ([4.10)1 ] . The moves of a tagged particle 
can be of two types, denoted (a) and (b). The sum of moves (a), where the particle 
follows a domain, is clearly asymptotically Gaussian, whereas moves of type (b), where 
the particle jumps from a domain to another, have a more complex distribution a priori. 
Numerical results point however toward a Gaussian dispersion in both cases. 

The crossover time between the coarsening regime and the equilibrium (or steady- 
state) regime scales as the equilibration time, which becomes very large at low temper- 
ature (£ ^> 1). It grows indeed as £ 5 for symmetric dynamics [see (|2.5|) ]. and as £ 4 for 
asymmetric dynamics [see (|2.6j) ]. 

Table El summarizes the low-temperature scaling behavior of the mean domain length 
and of the variance of the displacement of a tagged spin, in the steady state and in the 
nonequilibrium coarsening regime. For both types of dynamics, the dispersion of a tagged 
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Figure 8: Plot of the reduced third central moment (reduced cumulant) K% in the coarsening 
regime with asymmetric dynamics, against (t/£ 2 ) -1 . Symbols: data for several values of £. Full 
line: polynomial fit demonstrating the fall-off to zero. 



particle (measured by (r 2 ) c ) grows faster in the coarsening regime than at equilibrium. 
This behavior finds an explanation once interpreted in terms of a cage effect. In the 
coarsening regime, the effective cage felt by a spin consists of a single domain. Its size 
is therefore a typical domain length L, which grows as a power law, until it saturates 
to the equilibrium correlation length £. In the steady state, the situation is similar to 
that of the SEP or ASEP. The cage effect is a truly many-body effect, originating in 
the hard-core exclusion between particles in a high-density regime, and specific to the 
one-dimensional geometry. Remarkably enough, both the equilibrium behavior of (r 2 ) 
and the steady-state behavior of (r 2 ) c involve the ratio t/£. 



Symmetric 


t 


L 


(r 2 ) 


coarsening (Sec. 3.2) 
equilibrium (Sec. 3.1) 


e « t « e 


(t/a i/3 


(t/e? /3 


Asymmetric 


t 


L 


(r 2 ) c 


coarsening (Sec. 4.2) 
steady state (Sec. 4.1) 


e « t « e 
t»e 4 


(t/er 2 


t/e 
m 2/3 



Table 3: Dynamical regimes for low-temperature symmetric and asymmetric Kawasaki dy- 
namics, in the steady state, and in the nonequilibrium coarsening regime, with scaling behavior 
(powers of t and £) of the typical domain length, and of the displacement variance of a tagged 
spin. 
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The present work might be extended in several directions. Let us first mention the 
crossover between symmetric and asymmetric dynamics, in the presence of a weak bias. 
In the simple case of the mean domain length [3] , the essential ingredient of this crossover 
behavior is the gambler's ruin problem in the presence of a weak bias. Another facet of 
the present problem concerns the two-time correlation of the displacement, C(t + r, t) = 
{(r(t + t) — r(t)) 2 ). In the low-temperature coarsening regime, this correlation function 
is expected to exhibit a non-trivial dependence in the two-time plane. 
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